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P^h ■ Abstract 

A derivation of the Dirac equation in '3 + 1' dimensions is presented based on a 
master equation approach originally developed for the '1 + 1' problem by McKeon and 
,-Ch ' Ord. The method of derivation presented here suggests a mechanism by which the 

work of Knuth and Bahrenyi on causal sets may be extended to a derivation of the 
i -pH , Dirac equation in the context of an inference problem. 
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1 Introduction 

The Feynman Checkerboard (or Chessboard) problem [TJ [2] is a model from which the Dirac 
Equation[3] in '1 + 1' dimensions may be derived. Feynman's version of the problem was 
first published in his textbook on path integral methods[lJ. A combinatoric solution was 
published some years later [4J. That work, along with other combinatoric approaches, are 
critically reviewed and corrected elsewhere [5]. In addition to combinatoric methods of so- 
lution, Gersch[6] published a solution based on the correspondence between the Feynman 
Checkerboard model and the one-dimensional Ising model. This approach was developed by 
Ord and coworkers in a series of papers El ED]- In particular, the paper by McKeon and 
Ord [9] introduced a master equation approach for the solution of the '1 + 1' version of the 
problem which considered contributions from propagation events forwards and backwards 
in time. By imposing a causality constraint, the '1 + 1' Dirac equation emerged in a fairly 
straightforward fashion. The impressive variety of methods available for solving the Checker- 
board problem is complemented by the range of systems to which the basic model may be 
applied. As an example, work by Kholodenko|lll IT2] showed how the available approaches 
could be extended to, e.g., polymer dynamics and heterotic strings, arenas seemingly far 
removed from the original provenance of a textbook exercise. 

One shortcoming of the '1 + 1' Feynman Checkerboard model is that it does not account 
for spin. The history of spin within the context of a path integral approach is a checkered 
one, to coin a phrase. One possible approach was spearheaded by Schulman[T3]. In some 
sense, the success of the '1 + 1 model is due to the observation that spin seems not to exist 
as an independent concept in one spatial dimensionpj. Extensions to three dimensions have 
hitherto all seemed to founder on the difficulty of accommodating the required cr • p operator. 
The work presented here indicates one possible method for accomplishing this. 
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An intriguing application of this approach may be found in an extension of the investiga- 
tions of Knuth and Bahrenyi|!4] on causal sets, or posets, to a derivation of special relativity. 
The duality property of posets [15] suggested to the author that the forward and backward 
master equation approach exploited by McKeon and Ord [9] might offer some insights into 
how to extend the poset approach of Knuth and Bahrenyi|14] to a derivation of the '3 + 1' 
Dirac Equation. This is work in progress. The advantage of the master equation approach 
is that explicit expressions for real-valued quantities may be obtained at each time step, 
facilitating comparison to the poset approach. Complex amplitudes appear naturally as a 
result of a discrete (invertible) Fourier transform from time and space variables to a momen- 
tum space representation. At each stage, invertible, unitary transformations allow one to 
'follow the breadcrumbs' of a '3 + 1' Dirac equation back to the original, real- valued master 
equation. These assertions will be validated in the sequel. 



2 Master Equation Approach 

Define a transition rate a for changing direction on a line. It will be seen that this transition 
rate may be interpreted as a particle mass. The transition probability in time At will then 
depend on a At. McKeon and Ord [9] then write down a master equation for the probability 
amplitude of heading towards increasing values + or decreasing values — of x 

P ± (x,t + At) = (1 -aAt)P±(xT Ax,t) + aAtP T (x ± Ax,t). 

As McKeon and Ord showed[9j, the master equation for P±(x, t) may be iterated to obtain 

P + (x,t) = (! - a£t) n - R (aAt) R 
paths 

where n is the number of steps, and R is the number of reversals. Setting aAt = it and 
(1 — aAt) rs 1 reproduces Feynman's expression [1] for the probability amplitude. It is not 
clear from the analysis, however, that the substitution aAt = ie is a reasonable thing to do, 
although it does give the right answer. 

In order to address this issue, McKeon and Ord [9] developed a more versatile approach 
based on coupled master equations, allowing for the possibility of forward and backward prop- 
agation subject to a causality constraint. In a notation adapted from Ord and McKeon |10j 
one may write down coupled master equations for propagation forwards and backwards along 
the z axis 

Z±(z,t) = [l-(( + + (_)At]Z ± (z^Az,t-At) 

+( T AtZ±(zT Az,t + At) +(±AtZ T (z± Az,t - At) (1) 
Z T (z±Az,t + At) = [l-(( + + (_)At]Z T (z,t) + ( T AtZ T (z,t) + (±AtZ ± (z,t) (2) 
Z±(z,t) = Z T {z± Az,t + At), (3) 

where Z± is the forward time propagation probability towards larger z, i.e., Z + , or smaller 
z, i.e., Z_, Z is the backward time propagation probability, and Equation [3] is the causality 
constraint. In addition, (± is the transition rate for propagation towards larger (+) or smaller 
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(— ) values of z. McKeon and 0rd[9] note that iteration of coupled equations can be difficult. 
Instead, they inferred a differential equation from a short time expansion of Equations [1] and [2] 
subject to the constraint Equation [3] Performing a Taylor series expansion in Az = vAt 
and retaining terms only to order At, one finds the following differential equation for the 
difference of Z± and Z T 



±v 



'dZ± 


dZ T ~ 




~dZ± 


dZ T ~ 


dz 


dz 


+ 


. 9t 


dt 



+ [C+ + C-] [z± - z T ] = [c ± - Ct] [z? - z±] (4) 



Define A±(z,t) = exp([C + + £_] t) (Z±(z,t) — Z T (z,t)). One can interpret the integrating 
factor exp ([Ct + C-] t) as a chemical activity, which controls the 'concentration' of 'up' and 
'down' transitions. Substituting the definition of A± into Equation HJ one finds 



dA ± dA ± 
±u-^- + = (C± - Ct) A ± 



dz 



Of 



(5) 



In order to make further progress, it is useful to Fourier transform Equation to eliminate 
the z and t derivatives and work in the energy momentum representation. Set v = c and 
define Fourier amplitudes as follows 



A = exp(-i(pz - Et) /h)a±(p, E) 

p,E 

= 'y^ j ex-p(—i(pz + Et)/h)a±(p, E) 



(6) 
(7) 



p,E 



Substituting Equations [6] and [7] into Equation [5] and noting that if the result is to hold for 
all times over all values of z, then the coefficients a±(p,E) and a±(p,E) must satisfy the 
following constraints 



=F icpa± + iEa± 
^icpa± — iEa± 



h (C± - Ct) a T 
ft (Ct - Ct) «t 



(9) 



One may allow p and E to depend on z and t to account for space and time varying potentials, 
as will be shown. To be consistent in the order of At retained, it is sufficient to retain only 
the first term in derivatives of exp(— i(pz ± Et)) with respect to z and t, i.e., terms in dp/dz, 
dp/dt, dE/dz and dE/dt and higher are dropped in the limit At — > 0. 



3 Interpreting the Equations 

The forward time amplitudes a±(p,E) and the time reversed amplitudes a T encode the 
symmetry of the poset used in Knuth and Bahrenyi|14j. Note that (C± — Ct) ^ as the same 
sign as v in Equation [5j One may therefore rewrite (C± — Ct) = ^ lUJ where u is a positive 
definite constant. One may obtain a dimensionally consistent equation by defining hoj = mc 2 , 
the Compton energy of the particle. Jacobson and Schulman[4] have an insightful physical 
interpretation of mass as being proportional to a rate of making path reversals. That concept 
seems to be extended here to a notion of inertia, where the probability of making a particular 
class of path reversal Ct or C- depends on which direction the particle is traveling in. 
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Multiplying Equations [8] and M through by % and setting c = 1, one may write down 
matrix equations for a± and a± as follows 



E 
E 



a + 
a_ 



+ ma y 
- [PC 



ma,, 



(10) 
fill 



It is a straightforward exercise to show that the Equations [10] and [TT] satisfy the relativistic 
dispersion relation E 2 = p 2 + m 2 in a system of units where c = 1. 
Using the transformation 



$ = $ 



-i 



-i</>/2 











it is possible to rewrite the equations for a± = a±e Tl37T / 4 and a± = a±e ±i37f//4 in a convenient 
matrix form 



P 

-m 





-m 

~P 





) 
) 

■p —m 

m p 











= E 

















(12) 



subject to the constraint E 2 = p 2 + m 2 . Equation [12] is an eigenvalue equation in the form 
Hty = E^. Focusing on H, one may apply the following symmetry transformation 



s- 1 



10 
10 
10 
1 



to put H into the following form 
H = 

Note that Equation [13] is in the form 
H = 



p 





—m 








-P 





— m 


— m 





-P 








—m 





P 



(13) 



P -M 
-M -P 



Observe that M is proportional to the unit matrix and is invariant under unitary transfor- 
mations. Define 



P -> P" 



p" z 
-p" z 
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cos 6/ "2 sin 0/2 
sin 6/2 -cos 6/2 



P'z Px 
P'x -P'z 



and consider a rotation about the y" axis of the form 

Ug = U e l 

such that tan# = p' x /p' z . In this new coordinate system 

p" A p' = 

Now perform a rotation around the z' axis of the form 
such that tamp = p y /p x . In this new coordinate system p' z = p z . Thus 

p" 4- p' A p = Pa: _ ^ 

[Px + iPy ~Pz 

which properly encodes the operator cr • p. The following orthogonal transformation (R T 
indicates the matrix transpose of R) 



e iip/2 



1 -1 
1 1 



allows one to write 
RHR T -- 

which is equivalent to Dirac's time-independent equation in momentum space 



M cr ■ p 
cr ■ p —M 



[a ■ p + f3m] if> = Eil> (14) 

where the a and matrices are in the Dirac represent at ion [3]. Note that the components 
of ip are linear combinations of the Fourier amplitudes defined in Equations [10] and [TTJ The 
Fourier amplitudes may be traced back, ultimately, to the master equations for the forward 
and backward transition probability amplitudes defined in Equations Q] and |2] as all of the 
transformations are invertible. Given that the various a quantities are constantly being 
rephased and formed into new linear combinations, it is clear that they must be probability 
amplitudes and not probabilities per se. The picture that emerges is that a particle follows 
a stochastic trajectory in time and position, subject to the causality constraint introduced 
above. Manipulations of the relevant equations are facilitated by working in the momentum- 
energy representation. 
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4 Incorporating a Potential 



Given the assumption that p and E may be functions of position and time, one may reinter- 
pret the Fourier coefficients as canonical momenta. In this way, one would incorporate an 
electromagnetic four-vector as follows 

cAJ 



p — > p — 



E -> E- eA°/c, 



where r £ {1, 2, 3}. 

Under the assumption that 



/' = \l (Px - eA x /c) 2 + (p y -eA y /c) 2 + (p z - eA z /c) 2 
one may use the same series of steps used to derive Equation [TU to show that 
rmrT _ \ M tr-(p-eA) 

riri ti — , . - 

cr • (p - eA) —M 

where c = 1. Note that, expressed as a four-gradient, the four-momentum transforms as a 
covariant (lowered index) four-vector. With this observation the Dirac equation becomes 

[a • (p — eA) + /3m] ip = (E — eA )ip (15) 
One can recover the standard Dirac equation in the space time representation by substituting 
i) = J2 exp(i(p ■ x - Et))V(p, t) 

in Equation [T5j Note the argument of the exponent is a scalar, so if it is valid in one frame 
of reference, then it is valid in all reference frames. Then for i()(x., t) to be valid over all space 
and time, one requires 

eVWOM) (16) 



[a ■ (-iV - eA) + f3m] ip(x, t) , 

Equation [16] completes the derivation of the time-dependent Dirac equation in the presence 
of a potential from the master equation approach of McKeon and Ord[9j. In Feynman slash 
notation, one has 

{f - 4- + m)ip = 0. 

The derivation given here is fairly simple and straightforward compared to other ap- 
proaches to the Dirac equation in '3 + 1' dimensions based on the Feynman checkerboard 
problem [2] and references therein]. This may be due to the observation that the only essen- 
tial property of spinors used here is that the quadratic expression p\ + p 2 + p\ = p 2 can be 
'bilinearized' into the form 

Pz Px ~iPv] = ( 2 + 2 + 2\ [ 1 

P X + iPy 'Pz \ y P *^ P V^ P *) I 1 

a result which has much more to do with analytic geometry than any notions of 'quantum 
strangeness'. In order to complete the program sketched in the abstract, it is necessary 
to show how the master equation may be inferred from the poset approach of Knuth and 
Bahrenyi[14j. As noted, this is work in progress. 



Pz Px - ip v 

Px + iPy ~Pz 
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